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Abstract — We  present  an  algorithm  for  extracting  3D  canonical 
scattering  features  observed  over  sparse,  bistatic  SAR  apertures. 
The  input  to  the  algorithm  is  a  collection  of  noisy  bistatic  mea¬ 
surements  which  are,  in  general,  collected  over  nonlinear  flight 
paths.  The  ontpnt  of  the  algorithm  is  a  set  of  canonical  scattering 
features  that  describe  the  3D  scene  geometry.  The  algorithm 
employs  a  pragmatic  approach  to  initializing  featnre  estimates 
by  first  forming  a  3D  reflectivity  reconstruction  nsing  sparsity- 
regularized  least  squares  methods.  Regions  of  high  energy  are 
detected  in  the  reconstructions  to  obtain  initial  feature  estimates. 
A  single  canonical  feature,  corresponding  to  a  geometric  shape 
primitive,  is  fit  to  each  region  via  nonlinear  optimization  of  fit 
error  between  the  complex  phase  history  data  and  parametric 
scattering  models  using  a  modification  of  the  CLEAN  method. 
Feature  extraction  results  are  presented  for  sparsely-sampled, 
nonlinear,  3D  bistatic  scattering  prediction  data  of  a  simple  scene. 

Index  Terms — synthetic  aperture  radar,  bistatic  scattering, 
featnre  extraction,  radar  target  recognition 

I.  Introduction 

We  propose  an  algorithm  for  extracting  attributed  3D  fea¬ 
tures  from  SAR  data  collected  over  sparse,  nonlinear  apertures 
that  have  frequency,  azimuth,  and  elevation  diversity.  Practical 
SAR  flight  paths  provide  only  a  sparse  sampling  of  a  3D  SAR 
aperture,  resulting  in  high  sidelobes  and  aliasing  if  traditional 
backprojection  or  Eourier  imaging  is  used.  Recent  application 
of  regularization  techniques  has  shown  promise  in  removing 
these  image  artifacts  to  obtain  high  quality  3D  images  with 
regions  of  compact  energy  [l]-[3].  We  propose  a  feature 
extraction  algorithm  that  uses  the  regularization  approach  to 
generate  such  sparse  3D  images  for  detection  of  scatterers  and 
initialization  of  feature  parameter  estimates.  The  algorithm 
processes  fully-polarimetric  phase  history  data  and  outputs 
a  list  of  scattering  primitives,  characterized  by  their  type, 
location,  orientation,  and  size.  We  present  the  algorithm  for  the 
general  case  of  bistatic  SAR  apertures,  but  it  is  also  applicable 
to  monostatic  SAR,  as  a  special  case  of  the  bistatic  geometry. 

In  Section  II  we  present  parametric  scattering  models  for 
canonical  features  of  interest.  In  Section  III  we  discuss  the 
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feature  extraction  problem  and  outline  an  algorithm  for  extract¬ 
ing  features  from  polarimetric  phase  history  data.  In  Section 
V  we  apply  the  algorithm  to  simulated  data  for  an  arbitrary, 
nonlinear,  sparse  bistatic  SAR  aperture. 


II.  Bistatic  Canonical  Scattering  Models 


We  assume  that  scattering  from  an  object  can  be  represented 
as  a  sum  of  responses  from  canonical  scattering  centers.  A 
set  of  3D  canonical  models  have  recently  been  developed  for 
bistatic  measurements  [4],  [5].  These  models  include  a  plate, 
dihedral,  trihedral,  top-hat,  cylinder,  and  sphere.  In  this  section 
we  briefly  review  the  models;  details  on  their  derivations  are 
found  in  [4],  [5]. 

Using  high-frequency  approximations  [6],  we  represent  the 
received  radar  signal  S'(fc,A;0;r)  as  a  sum  of  scatterers 
S'g(fc,  A;  ©gjEq).  The  sum  may  be  written  as 

Q 

Sik,  A;  0;  T)  =  ^P^,(A;0,)Mr,(fc,  A;0,)e^'=^«A;0,)  (i) 

q=l 


where  k  =  is  the  wavenumber  for  frequency  /  and  speed 
of  light  c.  We  use  the  shorthand  notation  A  =  (j)r,  Or) 

to  represent  the  scattered  signal  dependence  on  the  transmitter 
(subscript  t)  and  receiver  (subscript  r)  azimuth  ((/))  and  ele¬ 
vation  (9)  locations.  The  term  is  the  2x2  polarization 
scattering  matrix  of  the  gth  scatterer,  which  takes  the  form 


F/3g— odd(A,  0g) 


-1  0 
0  -1 


(2) 


for  odd-bo\mce,  scatterers  or 

-  C0s(Ct  +  Cr)  sin(Ct  +  Cr)] 
sin(Ct  +  Cr)  C0s(Ct  +  Cr)J 

for  even-bounce  scatterers,  where  Cr  ^6  the  object 

rotation  angles  about  the  transmitter  and  receiver  antenna  lines 
of  sight  (LOS).  The  subscript  E^  in  (1)  indicates  the  type 
of  scatterer  (plate,  dihedral,  etc.)  and  encodes  the  particular 
functional  form  of  the  amplitude  response  Mr^.  The  vector 
0g  contains  physical  parameters  corresponding  to  a  reflector 
of  type  Eg,  including  {x,y,z)  location,  roll  (7),  pitch  (6), 
yaw  {(p)  orientation,  and  length,  height,  and/or  radius  size 
parameters,  as  applicable.  Einally,  the  term  AR  in  (1)  is  the 


-F/3g^even(A,  0g)  — 
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bistatic  differential  range,  which  may  be  approximated  under 
a  far-held  assumption  as 

AR{A;Qg)  =  ARr{A;Qg)+  Xq(pos  (pt  cos  0(-|-  cos  (j)r  cos  0^) 
-I-  Uq  (sin  (j)t  cos  9t  +  sin  (j)^  cos  9^)  +  Zq  (sin  9t  -I- sin  9^  (4) 

where  {xq,yq,Zq)  is  the  center  location  of  the  scatterer  and 
ARr{A;  &q)  is  a  range  offset  term  due  to  the  distance  of  the 
scatterer  surface  from  the  scatterer  center  location. 

Table  I  lists  canonical  shape  bistatic  scattering  models  as 
functions  of  relative  aspect  angles  and  (</>(,, 0(,)  that 

encode  both  the  radar  antenna  aspect  angles  and  the  object’s 
roll,  pitch,  and  yaw  orientation.  (See  the  model  derivations  in 
[4],  [5]  for  details  on  computing  the  relative  aspect  angles). 
Table  I  also  lists  a  calibration  scale  factor  ^cai,  which  when 
substituted  in  place  of  the  scalar  amplitude  term  A  in  Mr, 
scales  the  peak  scattering  intensity  maxdS'r,  P)  to  equal  the 
peak  radar  cross  section  (RCS)  of  the  object. 

Substituting  the  model  terms  from  Table  I  into  (1)  we  obtain 
a  description  of  the  scattering  response  of  a  set  a  rehectors  as 
a  function  of  the  data  collection  parameters  (radar  antenna 
locations  and  frequency)  and  object  parameters  (shape  type, 
location,  size,  and  orientation).  We  assume  that  the  radar 
antenna  parameters  are  known  and  that  the  object  parameters 
are  unknown.  In  the  sequel,  we  consider  estimating  the  object 
parameters  to  obtain  a  physical  description  of  target  structure. 

III.  Feature  Extraction 

Using  the  models  in  Section  II,  we  consider  the  problem 
of  extracting  scattering  features  from  a  set  of  noisy  bistatic 
phase  history  measurements,  Y{k,A).  Under  the  assumption 
of  white  Gaussian  measurement  noise,  maximum  likelihood 
parameter  estimates  of  the  scattering  parameters  can  be  found 
by  solving  the  least  squares  (LS)  problem: 

min 

©,r 

Note  that  (5)  is  a  combined  detection  and  estimation  problem. 
The  detection  problem  is  to  determine  the  scattering  types 
r  =  [Ti, . . .  jEat],  and  the  estimation  problem  is  to  solve  for 
the  parameter  vector  ©  =  [©i, . . . ,  ©at]  that  characterizes  the 
location,  orientation,  and  sizes  of  the  scattering  features.  In 
addition,  there  is  a  model  order  selection  problem  of  choosing 
N  in  (5);  ideally  N  =  Q,  the  true  model  order  in  (1). 

For  the  models  in  Section  II,  (5)  is  a  high  order,  nonlinear, 
non-convex  optimization  problem  that  requires  accurate  initial¬ 
ization  for  iterative  numerical  solvers  to  converge  to  the  global 
minimum.  Previous  work  utilizes  energy  compactness  in  2D 
SAR  images  to  initialize  model  order  and  scatterer  location 
and  size  parameters  [7],  [8].  However,  3D  SAR  flight  paths  are 
often  sparse  apertures,  and  the  corresponding  Fourier-based 
or  backprojection-formed  images  contain  high  sidelobes  and 
aliasing,  hindering  feature  segmentation  in  the  image  domain. 
For  example.  Figures  2(a)  shows  the  phase  history  for  an 
arbitrary  sparse  bistatic  aperture,  and  Figure  2(b)  shows  the 


N 


Y(fc,A)-^S„(fc,A;  ©„,r„) 


(5) 


corresponding  Fourier  image.  From  the  image,  it  is  difficult  to 
discern  the  number  and  locations  of  the  scatterers. 

Recent  application  of  regularization  techniques  [1],  [2]  has 
demonstrated  effective  formation  of  high-quality  sparse  3D 
images  [3].  Therefore,  we  propose  a  pragmatic  3D  feature  ex¬ 
traction  algorithm  that  first  generates  sparse  3D  images  using 
^i-Regularized  Feast  Squares.  These  sparse  images  are  used 
to  initialize  model  order  and  scatterer  location,  orientation, 
and  size  parameter  estimates.  Parameter  estimates  are  refined 
for  one  scatterer  at  a  time  using  a  modification  of  CFFAN 
[9].  The  algorithm  output  is  a  list  of  detected  features  and  the 
corresponding  model  parameter  estimates.  Key  steps  of  the 
proposed  feature  extraction  algorithm  are  summarized  below. 

A.  3D  Sparse  Image  Formation 

Fet  the  underlying  scene  reflectivity  (sparse  image)  be  given 
by  s  and  let  fl)  be  a  mapping  of  the  reflectivity  to  the  observed 
phase  history  at  the  observed  transmit  and  receive  angles  A. 
Then,  the  -regularized  least  squares  problem 

s  =  argmin  ||Y  —  thsll?  +  A  ||s||,  (6) 

S 

may  be  solved  to  obtain  a  sparse  image  solution  s  with  sparsity 
enforced  by  the  scalar  constraint  A.  For  our  SAR  imaging 
problem,  the  mapping  fl)  is  given  by  HF^,  where  f  is  a 
discrete  Fourier  transform  matrix,  H  is  a  projection  of  the 
phase  history  onto  the  flight  path  samples,  and  superscript  H  is 
the  matrix  Hermitian  operator.  We  follow  [1],  [2]  to  efficiently 
compute  the  fi -regularized  sparse  image  s.  A  sparse  solution 
that  is  jointly  optimal  across  all  polarizations  is  preferred,  but 
for  computational  reasons  we  independently  compute  sparse 
images  for  each  polarization  channel  and  denote  the  solutions 
as  Sp,  pG  {HH,HV,VH,VV}.  We  form  a  composite  image 

Sc  =  [|shh|^  +  |shv|^  +  |svh|^  +  |svv|^]  ^  (7) 

from  the  images  corresponding  to  the  four  polarization  chan¬ 
nels.  An  example  sparse  image  is  shown  in  Figure  3.  We 
assume  each  region  of  high  energy  in  the  composite  image  Sc 
corresponds  to  one  scatterer  and  use  the  image  to  initialize  the 
scatterer  location,  orientation,  and  size  parameters. 

B.  Image  Segmentation 

We  extend  the  Peak  Region  Segmentation  algorithm  in  [10] 
to  segment  three-dimensional  SAR  images.  We  implement  a 
region  merging  method  that  accounts  for  noise  ripple  across 
extended  radar  targets  but  still  separates  closely-spaced  scatter¬ 
ers.  Details  are  found  in  [11].  Figure  4  shows  that  application 
of  this  automated  segmentation  algorithm  to  the  sparse  image 
in  Figure  3  results  in  seven  distinct  scattering  regions. 

C.  Parameter  Initialization 

We  assume  that  each  segmented  image  region  corresponds 
to  one  scatterer.  Therefore,  we  initialize  scatterer  location 
estimates  for  the  plate,  dihedral,  and  trihedral  shapes  at  the 
centroid  of  the  image  segment  of  interest.  For  the  cylinder,  top- 
hat,  and  sphere,  the  location  of  energy  in  the  image  is  offset 
from  the  scatterer  location  by  the  shape’s  non-zero  radius.  For 
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Shape 
Type  r 

Bounce 
Type  /3 

Amplitude  Response  Mr(fc,  A;  ©) 

Calibration 

Acal 

Range  Offset  Ai?r 

Plate 

% 

odd 

A.  sine  [ A:  ^  (sin  +  sin  0^)]  sinc[/c  ^  (sin  cos  6'^  +sin  cos  0^)] 

(pi-,  G  [ ~  1"  7  f-]  )  ^  [ “  T  ’  i"] 

LH 

0 

Dihedral 

even 

A  sinc[A:iT  (cos  6'^  —  cos  9'^^  ]  sinc[fc  ^  (sin  (/)J.  cos  6'^  +  sin  cos 

G  [O.f] 

(e'+e'\  n’  ^  f’ f  ] 

|cos(  *2  ’■j,  5] 

2LH 

0 

Trihedral 

m- 

odd 

2^^Asinc  [fc//  (cos0^— cos0^)]  |^sinc[/c7T  (cos(^^— cos0^— cos(^J— cos^JJ] 

+  sinc[fc7T (cos((^^  “1"  ~cos  ^t)  ]] 

|sm(*2  '■-I-”  tani(^Ye,,G[0,tan  Yv2)]||  *2  ’’  lY L  0]| 

|cos(*2  ’■-l-I  tanY^j.  e,.e[tan  Yy2)’5]||  sin(  ;j,0,.G[O.  J]  j 

0 

Cylinder 

odd 

COS  (p'^  sinc[fc  ^  (sin  (p'^  cos  9'^  +sin  <p'^  cos  0^)]  ,  (p'^,  (p'^^  [“I"  ’  5] 

r  cos^  ^  (cos  (p[  +COS  ep'^'j 

Top-hat 

even 

1  sm  (  *2  '■]  .  e^,e^e[o,  j] 

sine  [A;Tf(cos  —  cos  0^)1  x<  \  ,  ,\ 

1  cos(  *2  '■  j  .  0t.e,G[L  '] 

rcos[  S^)(cos.,+cosO 

Sphere 

o 

odd 

AvTr 

1 

r  (cos  +  cos  9'^ )  cos  ^  cos  ^ 

+  r  sin^^C^^^  (sin  O'f.  +sin  9'^) 

TABLE  I 

Parametric  bistatic  scattering  models  from  [4],  [5]. 


these  shapes,  we  test  multiple  scattering  locations  near  the 
image  segment  centroid  to  obtain  low-cost  initialization  points 
for  numerical  estimation  of  the  object  location  and  radius. 
Length  and  height  parameters  are  initialized  for  all  shapes 
by  computing  the  extent  of  the  segmented  image  regions. 
Scatterer  orientation  is  initially  estimated  by  computing  the 
rotations  between  the  image  coordinate  axes  and  the  principal 
scattering  directions  of  each  segmented  image  region.  For 
details  of  the  parameter  initialization  scheme,  see  [11]. 

D.  Feature  Shape  Classification 

The  image  segmentation  and  parameter  initialization 
schemes  provide  information  about  the  location,  spread,  and 
orientation  of  the  scattered  energy  in  the  scene,  but  they  do 
not  identify  the  canonical  shape  type  associated  with  each 
region.  We  implement  four  shape  discrimination  methods  that 
test  the  data  polarization  response,  model  fit,  predicted  radar 
cross  section,  and  predicted  object  size. 

1 )  Polarization:  We  use  the  polarization  response  from  the 
sparse  images  obtained  from  (6)  to  classify  observed  features 
as  odd-bounce  or  even-bounce  scatterers.  We  extend  the  po- 
larimetric  classification  scheme  in  [8]  to  define  odd-bounce 
and  even-bounce  scattering  basis  vectors  for  the  bistatic  case 
as 


0 

^odd  — 

0 

-^even  — 

COs(Ct  +  Cr) 

0 

.sin(Ct  +  Cr). 

for  LOS  rotation  angles  (jt  and  [11]. 

The  bounce  type  of  the  ith  voxel  in  the  SAR  image  is 
estimated  by  finding  the  least-squares  distance  between  the 


basis  vectors  B  in  (8)  and  the  polarimetric  scattering  response 
vector  Si=[sHHi,  |(sHVi  +SvhJ,  SyvJ^-  Formally, 


bounce^  =  arg  min 

Si-B'pSi 

r  1 

72 

0 

1 

72 

0 

O' 

1 

Bp 

4’, 

/3C{odd,even} 

1 

L72 

-1 

72 

0_ 

The  example  in  Figure  5  shows  polarization  bounce  type 
estimates  for  the  image  segments  in  Figure  4  that  correspond 
to  the  correct  bounce  types.  Thus,  we  use  choose  the  sub¬ 
set  of  canonical  shapes  with  the  estimated  bounce  type  as  the 
candidate  scatterer  types  for  each  image  segment. 


2)  Model  Fit:  The  six  scattering  models  in  Table  I  have 
similar  functional  forms,  and  confusion  between  models  may 
arise  when  estimating  the  best  model  fit  to  the  data.  Therefore, 
we  fit  to  the  data  the  j  =  1,...,  J  candidate  models  whose 
polarization  matches  that  determined  by  the  polarization  clas¬ 
sification  in  (9)  and  make  a  list  F  of  best-fit  models  Fj  as 


F  =  arg 


er,- 


min  er 

r-  ^ 

^  -7 


<  Tfit 


(10) 


where  er  is  a  scaled  version  of  the  non-coherent  cost  metric 

^  3 

in  (13)  defined  as 


er,  =  20  logio 


1  Y^Y-2  Y^Si 


SS  Sr 

Tj  1 J 


YHy 


(11) 


for  K  data  samples.  Additional  model  discrimination  tests  are 
applied  to  the  best-fit  model  estimates  to  determine  the  final 
feature  shape  classification  decision. 
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3)  RCS  Consistency:  The  RCS  consistency  check  proposed 
in  [12]  verifies  that  the  estimated  feature  size  corresponds  to 
the  observed  scattering  intensity.  We  define  the  RCS  consis¬ 
tency  test  as 


C  =  arg  ■ 


20  log 


10 


Acai(T,g,f,rj-) 


<  TrCS 


(12) 


where  ^(Tj)  is  the  estimated  model  amplitude  corresponding 
to  the  scalar  A  in  the  Mr  terms  of  Table  I,  Acai(T,  f ,  Tj)  is 
the  calibration  amplitude  factor  defined  in  Table  I,  evaluated 
using  the  estimated  object  sizes,  and  (7  is  a  list  of  feature  types 
Tj  satisfying  the  above  inequality. 

4)  Object  Size  Test:  We  also  check  that  the  estimated  object 
sizes  are  reasonable  compared  to  the  scene  size — for  example, 
the  results  in  Section  V  check  for  object  sizes  between  0.01m 
and  Im. 

5)  Combined  Model  Discrimination  Tests:  We  implement 
the  three  discriminators  of  model  fit,  RCS,  and  size  to  choose 
the  best  model  for  the  data.  In  practice,  multiple  or  no  shapes 
may  pass  all  of  the  discrimination  tests,  so  we  combine  the 
discrimination  tests  as  follows.  In  the  case  that  one  or  more 
models  pass  all  three  discrimination  tests,  we  choose  the  model 
with  the  lowest  cost  given  by  (11).  In  the  case  when  no 
candidate  models  pass  all  three  discrimination  tests,  we  choose 
the  lowest  cost  model  that  passes  the  object  size  test.  If  no 
models  pass  the  size  test,  then  we  choose  the  lowest  cost 
model. 


E.  Feature  Estimation  with  CLEAN 

The  CLEAN  [9]  strategy  iteratively  estimates  and  removes 
signals  from  the  data  set,  starting  with  the  largest  energy 
contributor  to  the  observations.  At  each  iteration  i,  we  estimate 
the  ith  shape  model’s  parameters  as 

0.,  =  arg inin{Y(*)"Y(*) -2 abs{Y(*)'s,(fc,  A;  0^, Ti)} 

+  Sf(fe,A;0,,r,)S,(fc,A;0,,r,)},  (13) 

for  each  candidate  type  T^,  choose  the  best  match  according 
to  the  shape  classification  tests,  and  subtract  the  corresponding 
signal  estimate  from  the  data  to  obtain  a  new  set  of  observa¬ 
tions 

yb+i)  =Y«-S,(A;0„f,),  (14) 

repeating  the  process  until  a  desired  number  of  signals  or  a  de¬ 
sired  level  of  energy  is  modeled.  For  point  scatterers,  the  SAR 
image  provides  both  a  location  and  amplitude  estimate  for  the 
largest  contributing  point  at  each  iteration.  For  the  canonical 
shape  models  of  interest,  not  all  parameters  are  estimable 
from  the  image,  but  the  sparse  3D  images  help  to  initialize 
the  search  region  in  high-dimensional  parameter  space.  Note 
that  in  (13)  we  have  approximated  the  LS  solution  with  a 
non-coherent  cost  function  that  contains  fewer  local  minima, 
allowing  for  less  accurate  initial  estimates  [11].  Nonetheless, 
we  constrain  the  optimization  of  (13)  such  that  the  location 
and  orientation  estimates  are  near  the  initialization  points  and 


the  size  and  amplitude  parameters  are  positive-valued  scalars. 
We  fit  all  candidate  models  with  the  bounce  type  determined 
from  the  polarimetric  scattering  response.  Then,  we  follow  the 
feature  discrimination  scheme  described  in  Section  III-D5.  We 
subtract  the  selected  scatterer  model  phase  history  from  the 
observations  and  repeat  the  process,  from  image  formation  to 
parameter  estimation,  until  a  desired  number  of  scatterers  or 
a  desired  energy  level  is  modeled.  A  summary  of  the  entire 
feature  extraction  process  is  given  in  the  next  section. 

We  note  that  estimation  errors  may  result  in  high-energy 
residual  scattering  in  the  image  domain,  in  which  case  the 
CLEAN  algorithm  will  attempt  to  fit  multiple  canonical  shapes 
to  the  same  region.  However,  since  the  scatterers  considered 
are  well-separated  and  nearly  uncorrelated,  we  do  not  re-select 
voxels  in  the  image  for  which  we  have  already  estimated 
a  response.  This  modification  of  CLEAN  is  computationally 
practical,  though  estimation  errors  are  still  present  in  the  phase 
history.  The  impact  of  the  errors  on  estimating  subsequent 
high-amplitude  scattering  signals  is  small,  but  low-amplitude 
scatterers  may  be  difficult  to  estimate  if  the  error  terms  become 
large  relative  to  the  scatterer  of  interest. 

IV.  Feature  Extraction  Algorithm  Summary 

Given  phase  history  collected  on  (or  interpolated  to)  a 
Cartesian  grid, 

1)  Compute  composite  £i -regularized  LS  image 

2)  Segment  largest-energy  region  in  image. 

3)  Classify  image  segment  polarization  response  as  odd  or 
even  bounce  scatterer  type. 

4)  Initialize  parameters  for  each  canonical  model  whose 
bounce  type  corresponds  to  that  estimated  in  the  region. 

5)  For  each  candidate  model,  refine  parameter  estimates 
by  numerically  optimizing  non-coherent  cost  function 
between  observations  and  candidate  model. 

6)  For  all  models  within  Tfit  of  the  best-fit  model,  compute 
RCS  consistency  check  in  (12). 

7)  Check  that  size  estimates  are  within  acceptable  bounds. 

8)  Choose  estimated  model  to  be  best-fit  model  that  passes 
the  size  and  RCS  consistency  checks  according  to  the 
hierarchy  given  in  Section  III-D5. 

9)  Subtract  estimated  model  from  observed  phase  history. 
Repeat  Steps  1-9  until  meeting  a  user-specified,  maximum 
number  of  iterations  or  desired  level  of  modeled  energy. 

V.  Numerical  Results 

We  apply  the  feature  extraction  algorithm  of  Section  III  to 
simulated  bistatic  SAR  data  of  a  scene  comprised  of  seven 
canonical  shapes.  The  scene  is  depicted  in  Figure  1.  We 
simulate  four  polarization  channels  of  SAR  data  using  the 
model  signals  in  Table  I.  In  this  experiment,  the  receiver 
traverses  the  arbitrary  “squiggle  path”  shown  in  Figure  2(a), 
and  the  transmitter  is  fixed  at  =  (0°,30°),  near  the 

center  of  the  receiver  path. 

Since  the  SAR  phase  history  data  lies  on  a  sparse  aperture 
in  three-dimensions,  the  corresponding  Fourier  image  (Fig¬ 
ure  2(b))  contains  high  sidelobes  and  aliasing  that  make  it 
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difficult  to  discern  the  scatterers  in  the  image.  However,  the 
sparse  image  shown  in  Figure  3  contains  distinct  scattering 
regions  corresponding  to  the  canonical  shapes  in  the  scene. 
The  modified  Peak  Region  Segmentation  algorithm  locates 
seven  scattering  regions  of  interest,  shown  in  Figure  4.  For 
each  image  segment  in  Figure  4,  we  estimate  the  polarization 
bounce  types  from  (9).  The  estimated  bounce  types  are  shown 
in  Figure  5  and  agree  with  the  true  shape  types.  Using  the 
image  segmentation  and  polarization  bounce  results,  we  obtain 
initial  parameter  estimates  for  the  candidate  model  types  in 
each  region.  We  then  fit  canonical  shape  models  to  the  phase 
history  data  using  the  CLEAN  method. 

Table  II  lists  the  estimated  size  and  amplitude  parameters 
for  each  candidate  shape  along  with  the  results  of  the  model 
discrimination  tests  described  in  Section  III-D.  Shapes  selected 
by  combining  the  model  discrimination  tests  are  highlighted 
in  green,  indicating  correct  classification.  The  estimated  pa¬ 
rameters  for  the  selected  shapes  are  listed  in  Table  III,  along 
with  the  true  feature  types  and  parameters.  The  final  shape 
estimates  are  depicted  by  iconic  representations  in  Figure  6,  in 
which  the  icon  represents  the  estimated  feature  type,  location, 
orientation,  and  size.  In  the  example  shown,  all  seven  shapes 
are  correctly  classified.  The  parameter  estimates  for  each  shape 
converge  near  the  true  parameter  values,  except  for  the  dihedral 
height  in  image  segment  ^1,  which  is  estimated  to  be  more 
than  twice  the  true  height  of  the  dihedral.  The  estimated 
dihedral  does  fail  the  RCS  consistency  test,  which  correctly 
indicates  a  poor  amplitude,  or  as  in  this  case,  a  poor  size 
estimate  due  to  convergence  to  a  local  minima.  However, 
Figure  6  shows  the  estimated  scene  geometry  is  in  general 
a  very  good  match  to  the  true  scene  in  Figure  1. 

VI.  Conclusion 


Fig.  2.  (a)  Composite  phase  history  for  a  receiver  “squiggle  path”  and  a  fixed 
transmitter  at  =  (0°,30°)  near  the  center  of  the  receiver  aperture 

and  (b)  the  corresponding  HH  polarization  Fourier  image. 


y  (meters) 


Fig.  3.  Composite  sparse  image  for  spai'sity  constraint  A  =  25. 
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Fig.  4.  Image  regions  detected  in  Figure  3  by  the  image  segmentation 
algorithm  described  in  Section  III-B. 


We  have  presented  an  algorithm  for  extraction  of  3D  geo¬ 
metrical  features  from  sparse,  non-linear  SAR  apertures.  The 
feature  estimation  approach  is  an  extension  of  the  CLEAN 
algorithm  that  iterates  estimation  and  removal  of  canonical 
shape  signals.  We  use  regularization  to  form  sparse  3D  images 
that  allow  us  to  discern  scatterers  in  the  scene  and  initialize 
3D  canonical  shape  model  estimates.  We  have  shown  example 
results  for  the  proposed  feature  extraction  algorithm  for  an 
arbitrary  bistatic  SAR  aperture;  however,  the  algorithm  may 
also  be  applied  to  monostatic  apertures  (see  [11]  for  exam¬ 
ples).  Euture  work  will  apply  the  proposed  feature  extraction 
algorithm  to  SAR  data  of  vehicle  targets. 


y  (meters) 


Fig.  5.  Bounce  types  estimated  from  polarimetric  responses  of  segmented 
image  regions  shown  in  Figure  4.  Red  squares  represent  even-bounce  voxels 
and  green  triangles  indicate  odd-bounce  voxels. 


Fig.  1.  Scene  comprised  of  seven  canonical  shape  scatterers. 


Fig.  6.  Iconic  representations  of  the  estimated  features. 
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Image 

Segment 

# 

Candidate 

Shape 

Type 

Estimated  Size  and  Amplitude 

Discrimination  Test  Values 

Does  Shape  Pass  Tests? 

L 

H 

f 

Acal 

A 

201ogio(%l) 

NC  Cost 
dB 

Fit 

(Tfit  ^  5) 

Size 

[0.01,  l]m 

RCS 

('ntcs=5) 

1 

dihedral 

0.604 

0.694 

— 

0.839 

0.364 

7.244 

-68.808 

yes 

yes 

no 

1 

top-hat 

— 

0.000 

0.000 

0.000 

0.011 

—  oo 

-65.488 

yes 

no 

no 

2 

dihedral 

0.600 

0.301 

— 

0.361 

0.360 

0.036 

-68.599 

yes 

yes 

yes 

2 

top-hat 

— 

0.820 

0.253 

0.980 

0.246 

12.007 

-65.709 

yes 

yes 

no 

3 

plate 

0.600 

0.300 

— 

0.180 

0.180 

0.002 

-68.859 

yes 

yes 

yes 

3 

cylinder 

0.700 

— 

0.071 

0.186 

0.447 

-7.596 

-66.382 

yes 

yes 

no 

3 

trihedral 

— 

0.467 

— 

0.756 

0.117 

16.245 

-66.417 

yes 

yes 

no 

3 

sphere 

— 

— 

0.346 

0.346 

0.453 

-2.325 

-65.599 

yes 

yes 

yes 

4 

dihedral 

0.026 

0.532 

— 

0.027 

0.047 

-4.684 

-69.024 

yes 

yes 

yes 

4 

top-hat 

— 

0.499 

0.100 

0.376 

0.376 

-0.005 

-69.235 

yes 

yes 

yes 

5 

plate 

0.172 

0.000 

— 

0.000 

0.020 

—  oo 

-67.803 

yes 

no 

no 

5 

cylinder 

0.171 

— 

0.103 

0.055 

0.180 

-10.293 

-68.410 

yes 

yes 

no 

5 

trihedral 

— 

0.120 

— 

0.050 

0.050 

0.054 

-68.351 

yes 

yes 

yes 

5 

sphere 

— 

— 

0.521 

0.521 

0.617 

-1.459 

-66.691 

yes 

yes 

yes 

6 

plate 

0.117 

0.037 

— 

0.004 

0.013 

-9.558 

-67.050 

yes 

yes 

no 

6 

cylinder 

0.122 

_ 

0.135 

0.045 

0.104 

-7.367 

-67.041 

yes 

yes 

no 

6 

trihedral 

— 

0.088 

— 

0.027 

0.042 

-3.823 

-66.933 

yes 

yes 

yes 

6 

sphere 

— 

— 

0.750 

0.750 

0.750 

-0.002 

-68.490 

yes 

yes 

yes 

7 

plate 

0.501 

0.020 

— 

0.010 

0.024 

-7.581 

-67.598 

yes 

yes 

no 

7 

cylinder 

0.500 

— 

0.151 

0.195 

0.194 

0.039 

-67.618 

yes 

yes 

yes 

7 

trihedral 

— 

0.033 

— 

0.004 

0.008 

-6.477 

-65.699 

yes 

yes 

no 

7 

sphere 

— 

— 

0.000 

0.000 

0.147 

—  oo 

-65.656 

yes 

no 

no 

TABLE  II 

Estimated  shape  sizes  and  amplitudes  and  model  discrimination  results.  Highlighted  rows  indicate  the  shape  selected  for  each 

REGION  USING  THE  MODEL  HT,  OBJECT  SIZE,  AND  RCS  CONSISTENCY  DISCRIMINATION  TESTS. 


Segment 

Estimate 

Detected 

Parameters 

# 

Type 

Shape 

X 

y 

z 

7 

e 

L 

H 

r 

A 

1 

Truth 

dihedral 

0.25 

0.5 

0 

0 

0 

0 

0.6 

0.3 

— 

0.36 

1 

Final 

dihedral 

0.250 

0.499 

0.000 

0.066 

-0.317 

0.033 

0.604 

0.694 

— 

0.364 

2 

Truth 

dihedral 

0 

1.1 

0.25 

-10 

0 

0 

0.6 

0.3 

0.36 

2 

Final 

dihedral 

0.000 

1.100 

0.250 

-10.004 

0.000 

-0.002 

0.600 

0.301 

— 

0.360 

3 

Truth 

plate 

-0.1 

0 

0.3 

0 

-30 

0 

0.6 

0.3 

— 

0.18 

3 

Final 

plate 

-0.100 

0.000 

0.300 

0.133 

-29.999 

-0.000 

0.600 

0.300 

— 

0.180 

4 

Truth 

top-hat 

0.25 

-1.5 

0 

0 

0 

0 

— 

0.5 

0.1 

0.38 

4 

Final 

top-hat 

0.250 

-1.500 

-0.000 

-0.005 

-0.002 

0.000 

_ 

0.499 

0.100 

0.376 

5 

Truth 

trihedral 

0.75 

0 

0.1 

0 

0 

0 

_ 

0.12 

— 

0.05 

5 

Final 

trihedral 

0.748 

0.000 

0.106 

-1.222 

-0.375 

-0.768 

_ 

0.120 

— 

0.050 

6 

Truth 

sphere 

-1.5 

-1 

0 

— 

— 

— 

— 

0.75 

0.75 

6 

Final 

sphere 

-1.500 

-1.000 

0.000 

— 

— 

_ 

— 

— 

0.750 

0.750 

7 

Truth 

cylinder 

0.5 

2 

0 

0 

0 

-5 

0.5 

— 

0.15 

0.19 

7 

Final 

cylinder 

0.499 

2.000 

-0.001 

0.001 

0.912 

-4.999 

0.500 

— 

0.151 

0.194 

TABLE  III 

Canonical  shape  parameter  truth  and  final  parameter  estimates  for  the  detected  shapes  highlighted  in  Table  II. 
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